Quantum Interference between Impurities: Creating Novel Many-Body States in 

.s-wave Superconductors 
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We demonstrate that quantum interference of electronic waves that are scattered by multiple 
magnetic impurities in an s-wave superconductor gives rise to novel bound states. We predict that 
by varying the inter-impurity distance or the relative angle between the impurity spins, the states' 
quantum numbers, as well as their distinct frequency and spatial dependencies, can be altered. 
Finally, we show that the superconductor can be driven through multiple local crossovers in which 
its spin polarization, (s 2 ), changes between (s z ) = 0, 1/2 and 1. 
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Over the last two years, several beautiful experiments 
have studied quantum interference of electronic waves 
that are scattered by multiple impurities pi 0, S ^] . In a 
groundbreaking experiment, Manoharan et al. |tfl] used a 
corral of magnetic impurities on the surface of a metal- 
lic host to demonstrate that quantum interference can 
lead to the focusing of electronic waves into a quantum 
image. Moreover, using scanning tunneling spectroscopy 
(STS), Derro et al. were the first to observe four res- 
onance states in the local density of states (DOS) of the 
one-dimensional chains in YBa2Cu30g+x. These states 
were interpreted as arising from quantum interference of 
electronic waves scattered by two magnetic impurities |j) . 
Quantum interference effects were also studied in optical 
quantum corrals by Chicanne et al. and between im- 
purities located on quantum dots by Hollcitner et al. [|J. 
Some first theoretical work 0] has focused on impurity 
geometries in metallic systems similar to the one studied 
by Manoharan et al. In contrast, quantum interference 
in strongly correlated electron systems, such as super- 
conductors (with the exception of NbSe 2 R ) , charge- and 
spin-density- wave systems, or even semi-conductors, have 
not yet been addressed. However, the study of interfer- 
ence effects in theses systems involving spin impurities 
is not only of great fundamental interest, but might also 
possess important applications in the field of spin elec- 
tronics [|| and quantum information technology [El. 

In order to describe the properties of complex impurity 
structures such as quantum corrals, it is first necessary to 
understand interference effects associated with the pres- 
ence of few impurities. In this Letter we therefore con- 
sider two impurities embedded in a general s-wave super- 
conductor (SC). The presence of two magnetic impurities 
allows for a coupling of the bound states associated with 
a single impurity JToj , and gives rise to the emergence 
of novel many-body states. We show that the nature of 
these novel states, i.e., their quantum numbers, can be al- 
tered by varying the distance between the two impurities, 
Ar, or the relative angle between the directions of their 
spin moments, a. Moreover, we demonstrate that these 
changes are accompanied by local crossovers in which 



the spin polarization of the superconductor changes be- 
tween (s z ) — 0,1/2, and 1. We predict that the inter- 
play between the states' quantum numbers and the inter- 
impurity distance determines the distinct frequency and 
spatial dependence of the two-impurity bound states. Fi- 
nally, we discuss the implications of our work for systems 
with a larger number of impurities. 

Starting point for our calculations is the T-matrix for- 
malism |l0| which we generalized to treat the case of 
N impurities of spin S with non-magnetic and magnetic 
scattering potentials & In the following, we focus 
on the case N = 2, and, following Ref.[|o|, treat the 
impurity spins as classical, static variables, correspond- 
ing to the limit /3q — JS/2 = const, and S — > oo. In 
a fully gaped s-wave SC, this approximation is well jus- 
tified since no Kondo-cffcct occurs for sufficiently small 
coupling between the impurities and the delocalized elec- 
trons. Within this approach, any interaction between the 
impurities is only important to the extent that it deter- 
mines the angle, a, between the direction of the impu- 
rity spins. Within the Nambu-formalism and in Matsub- 
ara frequency space the electronic Greens function in the 
presence of N impurities is given by 

G(r,r',uj n ) = G (r,r',uj n ) 

N 

+ 2J G ! o(r,ri,a; n )T(ri ) r J -,w n )Go(r J -,r / ,w„) , (1) 

where the T-matrix is obtained from the Bethe-Salpeter 
equation 

T(n,rj,uj n ) = V ri 8 n>rj 

N 

+V ri ^G Q (r i ,ri,uj n )f(ri,r j ,u! n ) . (2) 
1=1 

In the case of two impurities 

V Tl = ^ {U\cfq + J1S0-3) T3 ; 

V ra = 2 + cos a + J 2 Sai sin a) t 3 ; 

Go 1 (k, iu n ) = [iLU n TQ - 6kT 3 ] CTQ + AkT 2 CT 2 . (3) 
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FIG. 1: DOS for the clean case (solid line), a single impurity 
(dashed line), and two impurities (dotted line) for /3o = 300 
meV (the DOS is shown at the impurity site). 



Here, V ri ,r 2 is the scattering matrix for the impurities 
located at ri and r 2 , respectively. Without loss of gen- 
erality, we take the spin of impurity 1 to be parallel to 
the i-axis, while that of impurity 2 is rotated from the 
z-axis into the z:r-plane by an angle a. Go(k, iu> n ) is the 
Greens function of the unperturbed (clean) system in mo- 
mentum space, and <Tj, Ti are the Pauli-matrices in spin 
and Nambu-space, respectively. Ui and Ji are the poten- 
tial and magnetic scattering strengths of the impurities. 
We consider a two-dimensional (2D) electronic system 
whose normal state dispersion is given by £k = k 2 /2m — ^ 
(h = 1), where ji = kp/2m is the chemical potential, and 
kp = 7r/2 is the Fermi wave- vector (we set the lattice con- 
stant ao = 1). The results and conclusions presented be- 
low are qualitatively robust against changes in the form of 
£k, the dimensionality of the s-wave SC, or the size of the 
momentum-independent SC gap, Ao- For defmiteness we 
set /i = 370 meV and m~ 1 /Ao = 15, but quantitatively 
similar results are obtained for m _1 /Ao = 30. The DOS, 
N(r,ui), presented below is obtained from a numerical 
computation of Eqs.(Q)-@ with N(r,u) = An + A 22 
and An (r, u) = —Ira Gu (r, oj + iS)/n. 

For a single magnetic impurity in an s-wave SC, the 

* .(12). 

T-matrix possesses poles at frequencies uj re 's , reflecting 
the presence of two bound states. The spectroscopic ev- 
idence for these bound states are two peaks in the DOS, 
as shown in Fig. [I], where we present the DOS obtained 
from Eqs. ([!])- (Q) at the impurity site; for comparison, we 
also plot the DOS of the clean system. These results 
are in general agreement with those of STS experiments 
p2[ , which provides further support for the validity of 
the T-matrix approach. Assuming for defmiteness that 
the impurity spin S||z and J > 0, we find that the bound 
state at luHI < (wres > 0), which we denoted by \p, 1) 
(\h, I)), is particle-like (hole- like) with spin along the — z- 
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r [a ] along R=(r,0) 

FIG. 2: (a) Spatial dependence of peak (1) and (2) (see Fig. [j]) 
in the DOS along the i:-axis (R = (r, 0)) for two magnetic 
impurities located at r = and r — 2 with parallel spins. 
Inset: Spatial dependence of the DOS at u = ±Ao- (b) Same 
as (a) but with impurities located at r — and r = 4. 



direction (+z-direction). 

We next consider two magnetic impurities with parallel 
spins, Ui — and Ji — J. For Ar — oo, the two sets 
of bound states given by \p, and \h,1,i) (i = 1,2) 
are degenerate. However, for Ar < oo, the probability 
that an electron scattered by one of the impurities is also 
scattered by the second one is non-zero. Hence, quantum 
interference of electronic waves that are scattered by both 
impurities leads to the formation of novel even and odd 
(or bonding and anti-bonding) states, \p, l) e ,o = (\p,l 
, 1) ± \p,l,2))/y/2, and similarly for the hole-like states. 
This picture is confirmed by the numerically computed 
DOS shown in Fig. [I] for two impurities located at ri = 
(0,0) and r 2 = (2,0), and /3 = 300 meV (the DOS 
shown is that on one of the impurity sites). As expected, 
the DOS exhibits four mid-gap peaks with peak (1), (2) 
corresponding to the particle-like states \p, |) e ,o an d peak 
(3), (4) to the hole-like states \h, T) e ,o- 

To determine which peaks in the DOS correspond to 
the even and odd states, we plot in Fig. ||a the spatial 
dependence of the particle-like states (1) and (2) along 
the x— axis with R = (r, 0) (the location of the impuri- 
ties at r = and r = 2 are indicated by arrows). Since 
the DOS of the odd states vanishes by symmetry at the 
midpoint between the two impurities, i.e., at r — 1, peak 
(2) and (1) correspond to the odd and even particle- like 
states, respectively. Note that their spatial dependence is 
remarkably different: while the odd state exhibits oscilla- 
tions well beyond the two impurity region, the even state 
is primarily confined to the region between the two impu- 
rities. This qualitative difference is associated with the 
(fci?r)-oscillations of the \p, |, z)-states. Since hp = tt/2 
and Ar = 2, the wave-functions of \p, 1,1) and \p, 1,2) 
are shifted by a phase Ac6 = fcfAr = jt outside the two- 
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FIG. 3: ui reB for the even (solid line) and odd state (dotted 
line) as a function of Ar for parallel impurity spins and (a) 
Pa = 300 meV; (b) (3 = 400 meV; and (c) [3 = 600 meV. 
The ranges of Ar with (s z ) = 1/2 are indicated by arrows, (d) 
Local spin polarizations, s z (r), along the f-axis, for a = 0. 



impurity region (r > 2 in Fig. ^a) and their spatial os- 
cillations are consequently out-of-phase. Thus, \p,i, 1) 
and \p, 1,2) interfere destructively for \p, |) e , and only 
weak spatial oscillations are observable in the DOS for 
r > 2. In contrast, \p, shows constructive interfer- 
ence of \p, (i = 1,2) and its spatial oscillations are 
enhanced for r > 2. By changing the inter- impurity dis- 
tance to Ar = 4 , with A<p = 2ir, the interference pattern 
between the even and odd states is exchanged, and \p, |) e 
(\Pil)o) now exhibits strong (weak) oscillations beyond 
the two impurity region, r > 4, as shown in Fig. ^b. 

In Fig. H we present the bound state energies of the 
states \p, |) e ,o as a function of Ar for a = (the frequen- 
cies of the \h, t) e ,o-states are obtained via uj res — > —uj res ). 
The amplitude of the oscillations in uj res decays as 
^ r -d/2 e Ar/£ j n d- dimensions, but since £ > vf/Ao ~ 20, 
the exponential decay is barely perceivable in Fig. |3|. The 
oscillations' spatial period, 1/kp, directly reflects that of 
the (fci?r)-oscillations in the wave-functions of \p, |) e , , in 
contrast, to the (2fci?r)-oscillations of the DOS. For cer- 



tain Ar, the splitting between the even and odd states 
vanishes and the two sets of bound states are decoupled. 
This decoupling arises whenever the bound state wave- 
functions of one impurity possess a node at the position 
of the other impurity. The same type of oscillations were 
obtained in Ref. for the case of NbSe2 . 

At small Ar, the bound state energy of \p,i) e crosses 
zero, and the state becomes hole- like, \h, \) e (at the 
same time, the state |/i,T)e transforms into \p, f) e ). We 
find that this zero-crossing of Lo res is accompanied by a 
crossover in the spin-polarization of the superconducting 
system which at T = is given by 



\Jd 2 rJ [Ah{t,uj) - A 22 {r 1 uo)} 



(4) 



This crossover is similar to the one predicted to occur 
when the scattering strength, /3o, of a single magnetic im- 
purity in an s-wave SC exceeds a critical value, (3 C [O, [ll| 
(for the band parameters chosen, we obtain (3 C « 460 
meV). At this point, the impurity breaks a Cooper-pair 
and forms a bound state with one of its electrons. Specif- 
ically, for S||z and J > 0, the spin polarization changes 
from (s z ) = for /3 < /3 C , to (a J) = 1/2 for /3 > f3 c . Sim- 
ilarly, for /3q > Pc/2, the system undergoes a crossover 
at Ar c , and for Ar < Ar c one electron of the broken-up 
Cooper-pair forms a single bound state with both impu- 
rities. For Ar — * 0, the DOS reduces to that of a single 
magnetic impurity with scattering strength 2(3q. Accord- 
ingly, \p, 1) moves towards the particle-hole continuum 
and vanishes for Ar = 0. Finally, a comparison of Fig. ^ 
and Fig. |||a, shows that, as expected, the spatially more 
confined bound state possesses a larger |w res | than the 
spatially more extended one. 

As (3q approaches (3 C from below, the number of 
crossovers increases, as shown in Fig. ||b for /3q = 
400 meV w 0.87/3 c . Due to the oscillatory behavior of 
Lu res , the bound state energy of the J. -state crosses zero for 
several value of Ar c . As a result, the spin-polarization os- 
cillates between (s z ) =0 and 1/2 and the electronic sys- 
tem can be tuned through multiple crossovers by varying 
Ar. The same tuning could also be achieved by keep- 
ing Ar constant and altering kp through changes in the 
doping level using recently developed field-effect transis- 
tor geometries Q . 

For (3q > @ c the superconductor exhibits a different 
crossover in which its spin polarization changes from 
(s z ) = 1 to (s z ) = 1/2. For Ar = oo, each impurity 
breaks one cooper pair and the spin polarization of the 
superconducting system is (s z ) — 1. As Ar decreases, 
one of the bound state energies crosses zero at least once, 
as shown in Fig. ||c where we present tu res for the |-states 
and 0o = 600 meV > (3 C . Accordingly, (s z ) changes from 
1 to 1/2. Since for Ar — > 0, the odd bound state vanishes 
by symmetry, the spin polarization reaches (s z ) = 1/2 for 
any value of /3q > /3 C . 
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FIG. function of the angle, a, between the direc- 

tion of the impurity spins for Ar = 1 and /3q = 400 meV. 



Changes in (s z ) are also reflected in the spatially re- 
solved spin polarization, s z (r) = J_ A duj(Aii — ^22), 
which due to the limited frequency integration is experi- 
mentally more easily accessible. In Fig. ^|d, we plot s z (r) 
along R = (r, 0) for two impurities located at r = 
and r = 2 and (s 2 ) = 0,1/2, and 1, corresponding to 
f3 = 300, 400 and 800 meV, respectively. For (s z ) =0, 
the spin polarization near the impurities is negative, as 
expected for S||z and J > 0. For (s z ) = 1/2, both im- 
purities form a single bound state with an electron in 
the I |) -state (see Figs. ||b). Thus, s z (r) is substan- 
tially increased at the impurity sites, but remains prac- 
tically unchanged at r = 1. In contrast, for (s z ) — 1, the 
electron from the second broken Cooper-pair joining the 
two-impurity bound state is in the | |) e -state, and con- 
sequently, s z (r) increases primarily around r = 1. Note, 
that for two impurities separated by Ar = 4, the first 
electron to join the two-impurity bound state is in the 
I J,) e -state, while the second one is in the | |) G -state, 
with corresponding changes in s z (r). 

The superconducting system can also be tuned through 
a crossover by changing the angle, a, between the two 
impurity spins, as shown in Fig. ||, where we plot cj res 
for all four bound states as a function of a (the impu- 
rities are located at ri = (u, 0) and T2 = (l,Uj). Since 



(indicated by the arrows on the right). This repulsion 
leads to the disappearance of all bound states for Ar — > 0. 

Finally, a non-zero U transfers spectral weight between 
the particle- and hole-like states and increases (3 C , but 
does not affect our above conclusions. Moreover, a self- 
consistent approach that allows for a gap suppression 
near the magnetic impurity does not change the quali- 
tative features of the DOS discussed above JO], in agree- 
ment with experiment 

The results presented above suggest that a supercon- 
ducting system with N impurities for which (3q > (3 C /N 
can be tuned through multiple crossovers with spin po- 
larizations ranging from (s z ) = to N/2, depending on 
/3o, the inter-impurity distances, and the angles between 
the spin moments. Work is currently under way to study 
these crossovers in more complex impurity geometries, 
such as quantum corrals, as well as the extensions to 
other host materials, such as unconventional SC, charge- 
density-wave systems, or semi-conductors j jl5| . 

In summary, we show that quantum interference of 
electronic waves scattered by two magnetic impurities in 
an s-wave SC gives rise to novel bound states. We predict 
that by varying the inter-impurity distance or the angle 
between the impurity spins, the states' quantum num- 
bers can be altered, and the SC can be driven through 
multiple local crossovers in which its spin polarization 
changes between (s z ) = 0, 1/2 and 1. 

We would like to thank A. Balatsky, J.C. Davis, A. 
de Lozanne, M. Randeria and A. Yazdani for stimulating 
discussions. 
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